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We analyze the contact and force networks in a dense confined packing of pen- 
tagonal particles simulated by means of the contact dynamics method. The 
particle shape effect is evidenced by comparing the data from pentagon pack- 
ing and from a packing with identical characteristics except for the circular 
shape of the particles. A surprising observation is that the pentagon packing 
develops a lower structural anisotropy than the disk packing. We show in this 
work that this weakness is compensated by a higher force anisotropy that 
leads to enhanced shear strength of the pentagon packing. With the polyg- 
onal shape of the particles, the strong force chains are mostly composed of 
edge-to-edge contacts with a marked zig-zag aspect. 



1 Introduction 

The microstructure of granular materials is generically anisotropic in two re- 
spects: 1) The contact normal directions are not random; 2) The force average 
as a function of contact normal direction is not uniform. The corresponding 
fabric and force anisotropics in shear are responsible for mechanical strength 
at the scale of the packing [11, 12, 13, 10]. The shear stress is fully trans- 
mitted via a "strong" contact network, materialized by force "chains" and 
the stability is ensured by the antagonist role of " weak" contacts which prop 
strong force chains [10, 14]. These features have, however, been for the most 
part investigated in the case of granular media composed of isometric (circular 
or spheric) particles. 

In this paper, we consider one of the simplest possible shapes, namely reg- 
ular pentagons. Among regular polygons, the pentagon has the lowest number 
of sides, corresponding to the least roundedness in this category, without the 
pathological space-filling properties of triangles and squares. We seek to iso- 
late the effect of edge-to-edge contacts on shear stress and force transmission 
by comparing the data with a packing of circular particles that, apart from 
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the particle shape, is fully identical (preparation, friction coefficients, particle 
size distribution) to the pentagon packing. Both packings are subjected to bi- 
axial compression simulated by means of the contact dynamics method. The 
presence of edge-to-edge contacts affects both quantitatively and qualitatively 
the microstructure and the overall behavior during shear. These contacts do 
not transmit torques, but they are able to accommodate force lines that are 
usually unsustainable in packings of disks. 



2 Numerical procedures 

The simulations were carried out by means of the contact dynamics (CD) 
method with pentagonal particles. The CD method is based on implicit time 
integration of the equations of motion and a nonsmooth formulation of mutual 
exclusion and dry friction between particles [15, 16, 17]. This method requires 
no elastic repulsive potential and no smoothing of the Coulomb friction law 
for the determination of forces. For this reason, the simulations can be per- 
formed with large time steps compared to molecular dynamics simulations. 
We used LMGC90 which is a multipurpose software developed in our labora- 
tory, capable of modeling a collection of deformable or undeformable particles 
of various shapes by different algorithms [18]. 




Fig. 1. Snapshots of a portion of the samples S2 (a) and SI (b) composed of circular 
and pentagonal particles, respectively. 



We generated two numerical samples. The first sample SI, is composed of 
14400 regular pentagons of three different diameters: 50% of diameter 2.5 cm, 
34% of diameter 3.75 cm and 16% of diameter 5 cm. The second sample S2, 
is composed of 10000 discs with the same polydispersity. Both samples were 
prepared according to the same protocol. The gravity was set to zero in order 
to avoid force gradients in the samples. The coefficient of friction was set to 
0.4 between grains and to with the walls. At equilibrium, both numerical 
samples were in isotropic stress state. The solid fraction was (f>o — 0.80 for SI 
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and 4>q = 0.82 for S2. The aspect ratio was h/l w 2, where h and / are the 
height and width of the sample, respectively. Figure 1 displays snapshots of 
the two packings at the end of isotropic compaction. 

The isotropic samples were subjected to vertical compression by down- 
ward displacement of the top wall at a constant velocity of 1 cm/s for a 
constant confining stress oo acting on the lateral walls. Since we are inter- 
ested in quasistatic behavior, the shear rate should be such that the kinetic 
energy supplied by shearing is negligible compared to the static pressure. This 
can be formulated in terms of an "inertia parameter" I defined by [19] 




where s = y/y is the strain rate, m is the total mass, and p is the average 
pressure. The quasistatic limit is characterized by the condition / <C 1. In our 
biaxial simulations, I was below 1CT 3 . 



3 Shear stress 

In this section, we compare the stress-strain behavior between the packings of 
polygons (sample SI) and disks (sample S2). For the calculation of the stress 
tensor, we consider the "tensorial moment" M l of each particle i defined by 
[20, 14]: 

<j = £/^, (2) 

cEi 

where /° is the a component of the force exerted on particle i at the contact 
c, r c p is the j3 component of the position vector of the same contact c, and 
the summation is runs over all contacts c of neighboring particles with the 
particle i (noted briefly by c € i). It can be shown that the tensorial moment of 
a collection of rigid particles is the sum of the tensorial moments of individual 
particles. The stress tensor <r for a packing of volume V is simply given by 
[20, 14]: 

- = ^E Mi = ^E/^ (3) 

iev cev 

where i c is the intercenter vector joining the centers of the two touching 
particles at the contact c. We extract the mean stress p = (<7i + a 2 )/2 7 and 
the stress deviator q = (a\ — u^jl, where a\ and ai are the principal stresses. 
The major principal direction during vertical compression is vertical, we then 
define the volumetric strain by : 

f V dV ' i (, AV \ 
£ p = / ITT = ln 1 + -7T < ( 4 ) 



where Vb is the initial volume and AV = V — Vq is the cumulative volume 
change. 
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Figure 2 shows the normalized shear stress q/p for the samples SI and S2 
as a function of shear strain e q = £1—62- For S2, we observe a hardening 
behavior followed by (slight) softening and a stress plateau corresponding to 
the residual state of soil mechanics [21]. For SI, we observe no marked stress 
peak. The residual stress is higher for polygons (~ 0.35) than for disks (~ 
0.28). The higher level of q/p for the polygon packing reflects the organization 
of the microstructure and the features of force transmission that we analyze 
in more detail below. 




0.1 0.2 0.3 0.4 



Fig. 2. Normalized shear stress q/p as a function of cumulative shear strain e q for 
the samples SI and S2. 



4 Fabric anisotropy 

The shear strength of dry granular materials is generally attributed to the 
buildup of an anisotropic structure during shear due to friction between the 
particles and as a result of steric effects depending on particle shape [22, 23, 
24]. 

The probability density function Pg(6), where 9 is the orientation of the 
contact normal n, provides the basic orientational statistical information 
about the granular the fabric. It is 7r-periodic in the absence of an intrin- 
sic polarity for n. Most lowest-order information is generally given by the 
second moment of Pe, called fabric tensor [25]: 

F a p = - f : n a (e)n fj (e)P e (6)d9 = ^ (5) 

where a and (3 design the components in a reference frame and N c is the total 
number of contacts in the control volume V. By definition, tr(F) = 1. The 
anisotropy of the contact network is given the difference between the principal 
values Fi and F2. We define the fabric anisotropy a by 
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a = 2(F X - F 2 ). (6) 

Figure 3 displays a polar representation of Pe(9) for the samples SI and 
S2 at e q = 0.3. We observe a nearly isotropic distribution for the pentagon 
packing in spite of shearing whereas the disk packing is markedly anisotropic. 
This is a surprising observation in view of the higher shear strength of the 
pentagon packing. 

The evolution of a is shown in Fig. 4 as a function of e q for SI and S2. 
The anisotropy stays quite weak in the pentagon packing whereas the disk 
packing is marked by a much larger anisotropy, increasing to ~ 0.3 and then 
relaxing to a slightly lower value in the residual state. The low anisotropy 
of the pentagon packing results from a particular organization of the force 
network in correlation with the orientations of edge-to-edge and vertex-to- 
edge contacts in the packing [26]. This leads us to consider the contributions 
of force and texture anisotropies to average shear stresses. 




Fig. 3. Polar representation of the Fig. 4. Evolution of the anisotropy a 

probability density function Pe of the with cumulative shear strain e q for the 

contact normal directions 8 for the samples SI and S2. 
samples SI and S2 in the residual state. 



5 Force anisotropy 

The angular distribution of contact forces in a granular packing can be rep- 
resented by the average force (/)(n) as a function of the contact normal 
direction n. We distinguish the average normal force (f n ){9) from the average 
tangential force (ft) (9). As P(9), these two functions can be represented by 
their Fourier expansions truncated beyond the second term [11, 12, 13, 10]: 



j (fn)(0) = (f){l + a n cos2(8-9 n )} 
\{ft)(0) = (f)a t sm2(9-9 t ) 



(7) 



6 



Emilien Azema, Farhang Radja'i, Robert Peyroux and Gilles Saussine 



where (/) is the average force, a n and a t represent the anisotropics of the 
normal and tangential forces, respectively, and n and t are their privileged 
directions. It is convenient to estimate the anisotropies from the following 
"force tensors" : 

' H ai =I(fn)(9)n a n p de, 

#S = J(ft)(6)n a tpde. 
o 

Then, we have tr(H^) = (/), and 

ir(n) „(n) 

a — 9 i ~ ^ 
a « - H {n> +H {n> ' 

at - z H (") +ff ("). 

where the subscripts 1 and 2 refer to the principal values of the tensors. 
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Fig. 5. Evolution of force anisotropies a„ (a) and at (b) as a function of cumulative 
shear strain e q in samples SI and S2. 




Fig. 6. Snapshots of normal forces in samples S2 (a) and SI (b). Line thickness is 
proportional to the normal force. 



Figure 5 shows the evolution of a n and a t with e q in samples SI and S2. We 
see that, in contrast to fabric anisotropies (Fig. 4), the force anisotropies in the 



Influence of particle shape on shear stress in granular media 7 

pentagon packing remain always above those in the disk packing. This means 
that the aptitude of the pentagon packing to develop large force anisotropy 
and strong force chains is not solely dependent on the global fabric anisotropy 
of the system. Indeed, due to the geometry of the pentagons, i.e. the absence of 
parallel sides, the strong force chains are mostly of zig-zag shape, as observed 
in Fig. 6, and the stability of such structures requires strong activation of 
tangential forces. This explains, in turn, the large value of at for pentagons, 
very close to a n , whereas in the disk packing at is nearly half of a n . 

The anisotropics a, a n and at are interesting descriptors of granular mi- 
crostructure and force transmission as they underlie the shear stress. Indeed, 
it can be shown that the general expression of the stress tensor Eq. (3) under 
some approximations leads to the following simple relation [13, 24]: 

- ~ -(a + a n + a t ), (10) 
P 2 

where the cross products a n a and a t a between the anisotropics have been 
neglected compared to the anisotropics, and it has been assumed that the 
stress tensor is coaxial with the fabric tensor Eq. (5) and the force tensors Eq. 
(8). Fig. 7 shows that Eq. (10) holds quite well both for pentagons and disks. 
This equation provides an amazingly good estimate of the shear stress from 
the anisotropies under monotonic shearing. 




0.1 0.2 0.3 0.4 

e 

q 

Fig. 7. Evolution of the normalized shear stress q/p for the samples SI and S2 with 
e q together with the corresponding predictions from its expression as a function of 
fabric and force anisotropies Eq. (10). 



A remarkable consequence of Eq. (10) is to reveal the distinct origins of 
shear stress in pentagon and disk packings. The fabric anisotropy provides a 
major contribution to shear stress in the disk packing (Fig. 4) whereas the 
force anisotropies are more important for shear stress in the pentagon packing 
(Fig. 5). In this way, in spite of the weak fabric anisotropy a, the larger force 
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anisotropies a n and a t allow the pentagon packing to reach higher levels of 
q/p compared to the disk packing. 



6 Conclusion 

The objective of this paper was to isolate the effect of particle shape on shear 
strength in granular media by comparing two similar packings with different 
particle shapes: pentagons vs. disks. We observed enhanced shear strength and 
force inhomogcncity in the pentagon packing. But, unexpectedly, the pentagon 
packing was found to develop a lower structural (fabric) anisotropy compared 
to the disk packing under shear. This low fabric anisotropy, however, does not 
prevent the pentagon packing from building up a strong force anisotropy that 
underlies its enhanced shear strength compared to the disk packing. 

This finding is interesting as it shows that the force anisotropy in a granular 
material depends on the particle shapes. This mechanism may be the predom- 
inant source of strength for " facetted" particles that can give rise edge-to-edge 
(in 2D) contacts allowing for strong force localization along such contacts in 
the packing. Since the fabric anisotropy is low in a pentagon packing, the role 
of force anisotropy and thus the local equilibrium structures are important 
with respect to its overall strength properties. With pentagon packings, we 
were able to demonstrate the nontrivial phenomenology resulting from the 
specific shape of particles as compared to a disk packing. We found a similar 
behaviour for other regular polygons (hexagons and higher number of sides) 
as well as polyhedral particles in three dimensions. 
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